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We study the possibility of fractional quantum Hall effects in HgTe quantum wells using exact 
diagonalization. Our results show that Laughlin states, the Moore-Read state, and the Read- 
Rezayi Z 3 state can all be supported. However, near the level crossing point (of the single-particle 
spectrum), the gap can be destroyed by Landau level mixing, and the Moore-Read state and the 
Read-Rezayi state dominate over their respective competing states only for wide wells. For smaller 
well widths the Moore-Read state crosses over to the composite fermion Fermi sea, while the Read- 
Rezayi state loses its dominance over the hierarchy state. 
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HgTe quantum wells have attracted a lot of atten¬ 
tions in the last several years. This is mainly due to 
the prediction!- and observation^ of the quantum spin 
Hall effect in this type of systems. However, it is nat¬ 
ural to ask whether it can exhibit other phenomena as¬ 
sociated with two-dimensional (2D) electronic systems. 
One of the most prominent such phenomena is the frac¬ 
tional quantum Hall effect. This phenomenon is interest¬ 
ing and important due to its strongly correlated nature 
and possible application in topological quantum compu¬ 
tation. Fractional quantum Hall effect in HgTe quantum 
well is interesting because of the strong spin-orbit inter¬ 
action in this material. The effect of strong spin-orbit 
interaction is two-fold. First, it affects the single-particle 
wavefunctions, so even though the electron-electron in¬ 
teraction is the Coulomb interaction just like in other 
2D materials, the matrix elements of the Coulomb in¬ 
teraction are different. This is equivalent to electrons 
with a parabolic band interacting with an effective in¬ 
teraction that is different from the Coulomb interaction. 
This is an important point because at some filling factors 
there are competitions between different incompressible 
states, even between incompressible states and compress¬ 
ible states, and which state is realized and how stable this 
state is depend on the details of the interaction. Second, 
the single-particle Landau level spectrum is also modified 
by the strong spin-orbit interaction. In particular, there 
are values of the magnetic field for which two Landau 
levels (LLs) cross each other. LL mixing could become 
important around such points. 


In this work we study the possibilities of different frac¬ 
tional quantum Hall effects in HgTe quantum well using 
exact diagonalization in the torus geometry^ . The ef¬ 
fect of spin-orbit interaction on the single-particle wave- 
function is taken into account with an eight-band cal¬ 
culation of the envelope function^. In the presence of a 
magnetic field, the envelope function has, with the axial 


approximation, the form 


Vw ,fc*(r) = exp(ifc x a;) 
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where k x is the single-particle momentum in the x di¬ 
rection [the difference in the prefactor from Ref. @ is 
due to a different choice of the vector potential, A = 
— Byx , with x being the unit vector in the 2 -direction], 
and e lkxX ip n are Landau wavefunctions for a parabolic 
band. By a parabolic band we mean a (strictly 2D) 

system with parabolic dispersion in the absence of a 

—2 

magnetic field, i.e. H pa ra = \p = (,Px,P y ) be¬ 
ing the 2D momentum], which becomes H pa ra ,mag = 
7^-(p+ - A) 2 in the presence of a magnetic field. The 
solution of the parabolic problem is standard^ and ip n = 
(— l) n [n£ 2 2 2n (n\) 2 ]~ie~^x~ ekx ' >2 H n (j — £k x ), where l is 
the magnetic length and H n is the Hermite polynomial. 
ip n appear in the ansartz m because they form a basis for 
the raising and lowering operators, in terms of which the 
p y operator can be expressed. [The phase (— l) n is chosen 
so that we have ap n = y/riip n -\ and a'ip n = y/n+ lp n +i 
with the correct sign, when the raising and lowering op¬ 
erators are defined as in Ref. 0, which in the current 
gauge have the explicit forms a = -±=(k x £— | — d y ), a) = 

-A=(k x £ — | + d y )] Because of the specific forms of the 
equations for the envelope functions in the axial approx¬ 
imation, with the choice of the indices of the ip n as in 
Eq. (JTJ), the y-dependence of each equation reduces to a 
common factor of <p n , leaving us a system of ordinary 
differential equations for the fi(z), which is then solved 
numerically. 

tp n with n < 0 are understood to be 0, so the lowest till 
is —2, and if -2 has only one nonzero component, contain- 
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ing ipQ. (The negative sign of the index in ip -2 is of no 
significance; one could have shifted the origin of till by 
defining = till+ 2 , then L would take nonnegative 
integer values. We are simply following the convention 
used in Ref. @.) Hence V ’-2 is similar to the lowest LL 
in systems with parabolic bands, in the sense that the 
only nonzero component of Eq. GD is proportional to ipo, 
except that here there is a finite width in the z direction, 
i.e. there is an additional factor fg°\z) (there is no such 
factor for strictly 2D electron gas). To avoid possible 
confusion, we repeat that the word “parabolic” refers to 
the dispersion in the absence of a magnetic field for a 
system whose eigenfunctions [e lkxX <p n {y)\ in the presence 
of a magnetic field we use as building blocks in Eq. ©■ 
This parabolic case is also the usual approximation for 2D 
electron gas. However, there is nothing parabolic in our 
current case otherwise. In particular, the confinement in 
the current case is due to the discontinuity of the con¬ 
duction and valence band edges, which are modelled as 
piecewise constant. These discussions about the single¬ 


particle wavefunctions are actually quite standard^. 

The other LL that is closest to the bulk gap is tpo, 
which has seven nonzero components, including three 
n = 0 components, three n = 1 components and one 
n = 2 component. Because of the n = 1 components 
it may have some similar behaviours to the n = 1 LL in 
conventional 2D electron gas, e.g. supporting the Moore- 
Read stated and the Read-Rezayi stated. 

If we ignore LL mixing, then the effect of the finite 
width and the mutli-component nature of the single¬ 
particle wavefunction can be absorbed into a form fac¬ 
tor in analogy to the four-band case in Ref. 13 Since 
we use a different geometry (torus instead of disk) and 
hence a different gauge (Landau instead of the symmet¬ 
ric gauge), it may be worthwhile to repeat some of the 
steps. The single-particle wavefunction is actually, after 
the magnetic translational invariance along the two edges 
of the system L\ = (Li, 0 ) and L 2 = (Z^z, A 21 /) (in the 
xy plane )& is taken into account, 


,j. = c£ 


kez 


_ 1 P 2 t> 2 

e 2i2 “ 


( 2 ) 


k x =2iv(j x -\-kNP)/L\ _ 


In the above expression, C is a normalization constant, M is as given in Eq. ©. N 4 , is the number of magnetic 

fluxes through the system, and j x = 1,2,replaces k x as one of the single-particle quantum numbers. The 
matrix element of the Fourier transform of the density is given by 


4 LLJi! ; n LL ,j xl (<2) = / 


( 3 ) 


Note that r and Q are three-dimensional vectors. In particular, Q = Q\\ + Q z z , where Q\\ is a vector in the xy plane 
and is limited to a lattice 
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( 5 ) 


After performing the 2D integral in the xy plane (over the parallelogram determined by L\ and L 2 ), we get (suppressing 
the label till) 

Uj x 2 ,j xl {Q) = |C | 2 J dzY^\M z )\ 2 eiQzZ LiS'j x 2 ^i + mexpl^^{ 2 j X 2 -m)h- L ni , ( 6 ) 

where the prime on the Kronecker delta means that the two subscripts are equal mod N#, L Ui are the 

Laguerre polynomials, and explicitly Q| = 47 r 2 77 -(~ m ^ x + n ) 2 + . Setting Q = 0 in the above expression 


gives us 


/1 00 
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Writing C = -^==C, we find that the normalization constant C is determined by 
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Defining 


F{Q hZ ) = \C\ 2 Y J \f^)\ 2 Ln i {^Q\( 2 ), 

i 

then the u matrix element can be written 


U jx2,jtcl ^jx2,jxl+m 
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The Coulomb matrix elements are, shortening the notations j x \ etc. to j\ etc., 
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( 9 ) 


( 10 ) 


where A is the system size (area in the xy plane), the prime on the summation sum indicates that the term m = h = 0 
should be omitted^, and e is the dialectric constant. Substituting Eq. © in Eq. ([ 101 ), we get 


Vj 1 , 32 ,33,ji ~ a 


p+oo 


d Q z 47re : 
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The Q z integral can be performed easily, 


r +oo 


dQ z ^QA^-Z!) e -Qn\zi-z 2 \ 


2tt Q 2 + Ql 


2Q || 


Then 


Vj 1 ,j 2 ,33,34 3ji +32,33 +34 


E ■ 

/c,n£Z 


27re 2 


eQ|| 


-7 r (Q||,«LL)e 17 ^ 


2—j2-33+kN$ d 


where the now means that h and j 2 — j$ + kN^ should not be simultaneously zero, and we have defined 

/ +oo 

dz 1 dz 2 e~ Q ^ Zl ~ Z2 ^F(Q^z 1 )F(Q^z 2 ). 

-OO 


( 12 ) 


(13) 


The effect of the complicated band structure of the single-particle problem has been absorbed into this form fac¬ 
tor, as advertised. The interaction part of the Hamiltonian are given in terms of the V elements as H ln t = 
\ 12j 1 j 2 j 3 j 4 Vj lt j 2 j 3 J 4 c ]i c j 2 c i 3 c i 4 ’ where the cl and c are creation and annihilation operators. 

In the vicinity of the level crossing point, we will take into account LL mixing by doing exact diagonalization with 
both LLs included. In this case we will need inter-LL matrix elements and a straightforward generalization of Ref. [l2] 
is carried out. Specifically, the inter-LL u matrix elements are 
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where we have defined 
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with the last factor being the associated Laguerre polynomial. The form factor for V)j LL1 j i;n:LL2 j 2; „ LL3 j 3; „ LL4i j 4 is 


jr(«LLl ,n-LL2 ,™LL3 ,™LL4) J 


/ +OO 

dz 1 dz 2 e- Q « lzi - Z2l F( nLL1 ’ nLL4 \-Qi h z 1 )F( nLL2 ’ nLL ^(Q lh z 2 ) 

-OO 


(17) 


[this replaces the F(Q\\,n ll) in Eq. (fl2l) . with the other factors unchanged]. Now the interaction part of the Hamil¬ 
tonian is -Hint = 2 5Zn LL l,J 1 ;n L L2,j2;nLL3,j3;nLL4,j4 ^"LLl ,ji;n L L 2 J2;nLL3.J3;nLL4 ,j4 C n LL1 C rtLL 2 J 2 C "LL3,j3 C "LL4 J4 ■ 


n 


Sample A, A^=15, 6 =—, n LL =-2 and 0 



FIG. 1. (color online) Gaps as functions of magnetic field for 
v = ±1/3, ±2/3 in sample A as calculated by exact diagonal- 
ization. Both till = —2 and till = 0 levels are included in 
the diagonalization. The magnetic flux through the unit cell 
is N<p = 15 flux quanta (by unit cell we mean the simulation 
cell, i.e. the parallelogram with Li and L 2 as two adjacent 
sides, not the unit cell of the atomic lattice). The aspect ra¬ 
tio, i.e. \L 2 /Li\ = 1 is chosen to be 1, and the angel between 
Li and L 2 (denoted 6 ) is chosen to be 7t/2, i.e. the unit cell 
is a square. 


We first consider sample A in Ref.0 and 1 /3-like states. 
For sample A, we do not know the electron density (which 
is due to a small residual doping), so we consider many 
different possibilities, which is equivalent to treating the 
density as a tunable parameter. Since the till = — 2 and 
till = 0 LLs cross each other at B c ~ 8.72T, LL mixing 
may be important. Therefore we include both LLs simul¬ 
taneously in the calculation^. Fig. [T] shows the gap as a 
function of the magnetic field for v = ±1/3, ±2/3 (y = 0 
corresponds to, at least in the noninteracting case, the 
lower of these two LLs completely filled, and the higher 
one completely empty). For B sufficiently away from B c 
on each side we find that the effect of LL mixing van¬ 
ishes. This is indicated by the fact that v = 1/3 and 
v = 2/3 have the same spectrum [i.e. particle-hole sym¬ 
metry within the till — —2 (0) level for B < (>)L? C ], 
which is also the same as that from a single-LL calcula¬ 
tion. Similarly, v = —1/3 and v = —2/3 have the same 
spectrum [particle-hole symmetry within the till = 0 
(—2) level for B < (>)L? C ]. In these regions we can safely 
use a single LL for the calculation and go to larger sys- 


Sample A, B=6.3T, v=l/3 or 2/3 (h ll =-2) 
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FIG. 2. (color online) Exact diagonalization spectra for sam¬ 
ple A of Ref. [i| at B = 6.3T for (a) v = 1/3 and 2/3 (the frac¬ 
tionally filled LL is tlll = —2), and (b) v = —1/3 and —2/3 
(the fractionally filled LL is till = 0). For all the data in this 
figure we used = 30 and |Li|/|L 2 | = 1, and each color 
corresponds to a different angle 6 € [^, £]. For the horizontal 
axis, k is the magnitude of k, which is a wavevectorlike many- 
body quantum number characterizing the relative motions of 
the electrons^ (and has no relation with the k x earlier in the 
discussion about the single-particle wavefunction). 


tern size (i.e. larger N ^). The spectra for B = 6.3T are 
shown in Fig. [2] as an example. The neutral gaps are de¬ 
termined by the minima of the magnetoroton branches 
(the low-lying curve in each subfigure). We note in pass¬ 
ing that the dispersion of the magnetoroton branch can 
be measured experimentally^. When B is closer to B c , 
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Sample A, v=2/3, B= 7.3T, A^=15, 6>=- 


CN I %- 
^ I UJ 

< 


0.15 


0.10 


0.05 


0.00 




FIG. 3. (color online) The spectrum for sample A with v = 
2/3, B = 7.3T, and N$ = 15. Both till = —2 and till = 0 
LLs are included in the diagonalization (the total number of 
eletrons in these two LLs is then 25). The dots are levels with 
(N- 2 ) ~ 10 electrons in the till = —2 LL (the higher level 
at this B). The dashes are levels with ( N -2 ) « 11 electrons 
in the till = —2 LL (i.e. one electron is promoted to this LL 
from the till = 0 LL, leaving a hole in the latter). 


for each filling factor there is a narrow region (on each 
side of B c ) where the gap is determined by inter-LL exci¬ 
tations (excitations where some electrons are excited to 
the higher LL from the otherwise filled lower LL) and 
is reduced compared to the magnetoroton minimum, see 
e.g. Fig. [3] For B even closer to B c the gap is totally 
destroyed by LL mixing. 
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Next we consider v = 1/2. To be concrete, we fix 
the electron density to 4.2 x 10 n cnC 2 (as for sample 
B in Ref. jfj). We consider several different well widths 
w. For the density we choose, v = 1/2 corresponds to 
B = 34.75T and is to the right of the level crossing point 
between till = — 2 and till = 0 for all the well widths 
and therefore the fractionally filled LL is the till = 0 
level (and for all well widths this level is sufficiently far 
away from other levels for LL mixing to vanish, so we will 
only keep this one level in the exact diagonalization). 

First we look at w = 8nm as for sample B of Ref. 
Unless otherwise specified the geometry of the unit cell 
is hexagonal. For N e = 10 , the ground state (GS) has 
k = 5<fi [<fi and cfc are the basis vectors defined in Eqs. (pT|) 
and ©] and the two fc’s related by rotational symme¬ 
try, i.e. 5q 2 and 5(<fi + q 2 ) [i.e. the GS is three-fold 
degenerate (not counting the trivial two-fold center-of- 
mass degeneracy£ii£)] . These three fc’s correspond to the 
centers of the edges of the Brillouin zone. This behavior 
is consistent with the Moore-Read Pfaffian (Pf) stated 
[or the anti-Pfaffia n 17 ' 18 (APf), the particle-hole conju¬ 
gate of the Pf]. The overlap with the Pfafhan state is 
considerable at 0.62 (the overlap with the APf is the 
same in the absence of LL mixing). In Ref. [l9| it was 
argued that since the Coulomb state is particle-hole sym- 


n 
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FIG. 4. (color online) Overlaps between the Coulomb GS 
in HgTe quantum well with the Pf/APf, the symPf and the 
CFFS for v = 1/2 and (a) N e = 10, (b) N e = 12, and (c) 
N e = 14, as functions of the well width w. 
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v= 1/2 (flLL=0), e=- 



v=3/5 (n LL =0), 0=-,lV e =15 
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FIG. 5. (color online) Gaps for i/ = 1/2 as a function of well 
width for N e = 10,12,14. 
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FIG. 6. (color online) Spectrum for v = 1/2, w = 32nm, with 
N e = 16. The blue dots are for the hexagonal unit cell, while 
the black dashes are for the square unit cell. The scheme of 
Ref. [Ill is used to fold the Brillouin zone to a quarter of the 
original size to account for the pairing nature of the Pf/APf 
state. In particular the GSs are mapped from the centers 
of the edges in the original Brillouin zone to the center of 
the folded one (k = 0). For the hexagon unit cell, the GS’s 
overlaps with the Pf/APf and the symPf are 0.50 and 0.86, 
respectively [the CFFS for this system size and geometry has 
k = — 451 — 4^2 (and fc’s related by rotational symmetry) in the 
unfolded Brillouin zone, which is different from the Coulomb 
state here, so there is no overlap]. For the square unit cell, 
the overlaps between the GS and the CFFS, the Pf/APf, and 
the symPf are 0.50, 0.55 and 0.86, respectively. 


metric in the absence of LL mixing, while the Pf state is 
not, we should consider the particle-hole symmetrized Pf 
(symPf) instead. In the current case, the overlap with 
the symPf is 0.81. However, this does not necessarily 
mean that the Coulomb state is the symPf. The other 
possibility is the composite fermion Fermi sea (CFFS)^, 
also known as the Rezayi-Read state (not to be confused 
with the Read-Rezayi state to be discussed below). There 


FIG. 7. (color online) The overlap between the GS and the 
H-CF state and the total overlap of the GS and lowest excited 
state (at k = 0) with the Z 3 doublet for v = 3/5 as functions 
of the well width. The inset is the splitting between the GS 
and the lowest excited state (at k = 0). 


is no known parent Hamiltonian, i.e. a Hamiltonian for 
which this state is exactly the GS, so to calculate the 
overlap with this state one would need to use Monte 
Carlo method, which would take too long for the system 
sizes we are using here. Fortunately, it is known (from 
Monte Carlo calculations with smaller system sizes) that 
the Coulomb state in the lowest LL in the purely 2D 
case [i.e. iL P ara,mag given below Eq. Q])] has overlap with 
the CFFS that is practically l 19 i 20 , so assuming this is 
true also for larger system size, one can approximate the 
CFFS with the lowest LL Coulomb stat o 21 ’ 22 (the latter 
has the added advantage of being exactly particle-hole 
symmetric, while for the former the particle-hole sym¬ 
metry is only approximate^). With this approximation, 
we find that in the current case the overlap with the 
CFFS is 0.99, much higher than that with the symPf. 
For N e = 12 the GS has k = 4<fi + 2 q2 (and some other 
values related by rotational symmetry) and is not con¬ 
sistent with the Pf/APf/symPf. The overlap with the 
CFFS is 0.99. 

For N e = 14 the GS has k consistent with the 
Pf/APf/symPf, the overlaps with the Pf/APf, the symPf 
and the CFFS are 0.39, 0.58 and 0.97, respectively. The 
results for all three system sizes suggest that at this value 
of well width (w = 8nm) the GS is the CFFS. How¬ 
ever, when we increase the well width, for N e = 10 
and 14, the overlaps with the Pf and the symPf in¬ 
creases while the overlap with the CFFS decreases, and 
for w > w c ~ 15nm, the overlap with the symPf is larger 
than those with the CFFS, see Figs. |4(a~)] and 

For N e = 12, the k of the GS changes from that cor¬ 
responding to the CFFS to that corresponding to the 
Pf/APf/symPf at around w = 21nm [Fig. |4(b)] . These 
behaviours suggest a crossover from the CFFS at small 
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FIG. 8. (color online) Exact diagonalization spectrum for 
v = 3/5 and w = 44nm with N e = 18 


well width to the symPf at large well width. In Fig. [5] 
we show the well width dependence (and system size de¬ 
pendence) of the gap. In Fig. [6] we show the spectrum 
for N e = 16 (both hexagon and square unit cells) at 
w = 32nm. The GS k values and the overlaps for this 
system size also suggest the symPf state. 

Finally, we consider v = 3/5. This filling factor is of 
particular interest because the GS could be in the uni¬ 
versality class of the Z 3 parafermion stated, which can 
support universal quantum computation^. An compet¬ 
ing state for this filling factor is the hierarchy stat o 24 ’ 25 , 
which can also be described as an integer quantum Hall 
effect of composite fermions^. As for v = 1/2, we fix 
the density to 4.2 x 10 11 cm -2 and vary the well width. 
For this density zx = 3/5 corresponds to B = 28.96T, 
for which the n ll = 0 LL is also far away from other 
LLs, so we only include this one LL in the exact diag¬ 
onalization. For w = 8nm, the overlap between the GS 


and the hierarchy-CF (H-CF) state is 0.88. The GS for 
the Z 3 state is two-fold degenerate (the GS k is equal to 
0). For HgTe quantum well the GS is nondegenerate (not 
including the trivial five-fold center-of-mass degeneracy) 
and the splitting between the Coulomb GS and the lowest 
excited state (at k = 0) is 0.043|y, which is quite large. 
However, if we ignore this splitting and calculate the total 
overlap of the GS and the first excited state with the Z 3 
doublet^, i.e. *|GS)| 1 2 3 4 5 6 * 8 9 +|(Z 3 ,f|FES)| 2 ), where 

\Z 3 , 1) and \Z 3 , 2) are orthonormal basis states of the Z 3 
GS subspace, |GS) and |FES) denote the GS and the first 
excited state (at k = 0) for the HgTe quantum well, we 
get 0.92. Hence for w = 8nm neither the H-CF state nor 
the Z 3 state is clearly favoured. As w is increased, the 
overlap between the GS and the H-CF state decreases, 
the splitting between the GS and the first excited state 
decreases, and the total overlap of these two states with 
the Z 3 GS doublet increases, see Fig. [TJ Therefore, at 
least for the larger well widths, the Z 3 state is a viable 
candidate for the GS at v = 3/5. In Fig. [H] we show the 
spectrum for w = 44nm and N e = 18. 

To conclude, our exact diagonalization results show 
that HgTe quantum wells are capable of supporting frac¬ 
tional quantum Hall effects. The filling factors we studied 
include v = ±1/3, ±2/3,1/2, and 3/5. For v = 1/2 we 
find a crossover from the composite fermion Fermi sea at 
small well width to the particle-hole symmetrized Moore- 
Read state. For v = 3/5, the Read-Rezayi Z 3 state dom¬ 
inates over the hierarchy state at large well width but 
gradually loses this dominance when the well width is 
decreased. We also find the effect of Landau level mixing 
can be important near the level crossing point and can 
totally destroy the gap. We note that our study does not 
take into account disorder, which can reduce the gaps 
significantly and make the experimental observation dif¬ 
ficult. 
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